Phase behavior of repulsive polymer-tethered colloids 
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We report molecular dynamics simulations of a system of repulsive, polymer-tethered colloidal 
particles. We use an explicit polymer model to explore how the length and the behavior of the 
polymer (ideal or self-avoiding) affect the ability of the particles to organize into ordered structures 
when the system is compressed to moderate volume fractions. We find a variety of different phases 
whose origin can be explained in terms of the configurational entropy of polymers and colloids. 
Finally, we discuss and compare our results to those obtained for similar systems using simplified 
coarse-grained polymer models, and set the limits of their applicability. 



I. INTRODUCTION 

Understanding how colloidal particles spontaneously 
organize into ordered macroscopic aggregates is a long- 
standing challenge that has recently acquired an extra 
degree of complexity. In fact, advances in particle syn- 
thesis [Il-IH have opened the way to the production of 
colloidal particles that are anisotropic both in shape and 
surface chemistry. This provides an unlimited number 
of building blocks that can spontaneously assemble into 
an unprecedented variety of structures with potentially 
novel functional, mechanical, and optical properties. 

The effect of the anisotropy of nanoparticles on their 
macroscopic ordering can be addressed in terms of (a) 
the form of the inter-particle interaction, and (b) their 
shape. Not surprisingly, for both cases there is ample 
evidence (see for example [6l-[l5j and references therein) 
of a strong correlation between the physical properties 
of the components and those of the resulting aggregates. 
This phenomenology must be thoroughly explored as it 
may lead the way to a rational design of the components 
to target desired macroscopic structures. 

Here we focus on the role of particle shape. Specifically, 
we study the phase behavior of a particularly interesting 
class of deformable particles that is obtained by graft- 
ing a single long chain to a each colloid. What makes 
this hybrid colloid intriguing is that, because of the flex- 
ibility of the polymer, the overall shape of the particle 
is not fixed, but can be spontaneously altered depending 
on the specific thermodynamic states imposed on the sys- 
tem. A few experimental realizations of particles compat- 
ible with our model have recently been synthesized |16l - 
[Tsl . The dual nature of these nanoparticles may open 
the door to exotic self-assembled structures that are not 
typically seen in systems of nanoparticles with intrinsic 
(invariable) shape, and its key elements are very similar 
to those of asymmetric diblock copolymers. 

Unlike recent experimental and theoretical studies on 
particles coated with dsDNA, which can form com- 
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plex networks between the particles via linker-mediated 
dsDNA-dsDNA interactions jl9l - [21| . no explicit attractive 
forces are introduced in our system. As a result, any 
phase described in this paper will be mostly driven by 
a nontrivial balance between the configurational entropy 
of the colloids and that of the chains. 

In this paper we use molecular dynamics simulations 
to understand the phase behavior of a system of repul- 
sive, polymer-tethered colloidal particles. Specifically, we 
consider a system in which each colloid is connected to 
one of the end groups of a single polymer, and we study 
how different structures emerge depending on the poly- 
mer length. Furthermore, we explicitly analyze both the 
case of ideal and self-avoiding polymers. 

Capone et al. 22| have recently analyzed the phase 
behavior of model of di-block copolymers. Their study 
bears similarities with our work, as some of the self- 
assembled structures are common to both systems. How- 
ever, the two systems differ in the way the polymers 
are modeled. Here we use an explicit bcads-and-springs 
model, while spheres with a soft potential were employed 
by Capone et al. Though computationally expensive, 
our model enables us to gain a detailed understanding 
of the mechanisms behind the nontrivial phase behavior 
emerging in this system, and this choice will turn out to 
be quite critical when considering the case of non-ideal 
polymers. To the best of our knowledge this is the first 
computational study that explicitly accounts for the in- 
ternal degrees of freedom of the polymer for this partic- 
ular system. 



II. MODEL 

We model the polymer-tethered colloids as a polymer 
of iV + 1 monomers, with N monomers of diameter cti, 
and monomer N + \ , representing the colloidal particle, 
of diameter (72 . See Fig. [T] for a depiction of the particle. 
In this model the iV*'' monomer is not constrained to be 
at a specific location on the surface of the colloid, but 
can freely diffuse on it; constraining this monomer would 
yield the same equilibrium properties. 

Excluded volume interactions between any two par- 
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FIG. 1: Schematic representation of our model for a hybrid, 
polymer-grafted colloid. The first monomers of diameter 
o"! represent the chain, whereas the A'^ + monomer of di- 
ameter (T2 represents the colloidal particle. Both particles are 
assumed to have equal mass. 

tides in the system are enforced via a purely repulsive 
shifted-truncatcd Lennard- Jones potential 
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The indices i,i G {1,2} indicate the identity of the par- 
ticle (polymer or colloid, respectively.) (Jij = {01+02) /'i- 
ri,j is the distance between the centers of mass of any two 
particles. Finally, e^j- = lOk^T V i, j when considering 
self-avoiding polymers, and we set ei^i = Ok^T for the 
case of ideal polymers. 

In each hybrid colloid, particles are linearly connected 
via the harmonic spring potential 

Ut+i = ks^^,^+l ~ rl^i)' (2) 

where kg = WOk^T is the spring constant, and r'^ = 
cr,^i+i + (Ti/2 is the equilibrium distance. 

We perform NPT molecular dynamics simulations us- 
ing the LAMMPS (Large-scale Atomic/Molecular Mas- 
sively Parallel Simulator) package [2^. Pressure and 
temperature are kept constant by means of a Nose- 
Hoover thermostat [23| and barostat [25| with additional 
drag terms, with coefficients = l''o~^ ^'^'-^ " ^'''0^ 
respectively (tq is the reduced time unit), to damp the 
dynamics and suppress large temperature and pressure 
oscillations. 

The simulation box is a cuboid with periodic bound- 
ary conditions and, for pressure control, we use decoupled 
box lengths in each of the three Cartesian coordinates. 
This allows box aspect ratios to vary slightly to accom- 
modate crystalline structures. 

The system initial configurations are prepared by per- 
forming NVT simulations in the gas phase. Once the 
system is equilibrated, and the initial pressure Pq is ex- 
tracted from the thermalized configurations, we slowly 
ramp the pressure to the desired value Pi starting from 
Pq (all pressures referred in this paper are rescaled with 
respect to the colloidal interaction energy £22 and the 
colloidal diameter 02)- Each subsequent simulation per- 
formed at a constant pressure P, starts from the ther- 
malized configuration at pressure P^-i (P^-i < Pi). The 



dimcnsionless AP was typically set to O.OI, and was re- 
fined or extended depending on the distance from the 
transition point. This procedure ensures that the chains 
have the time to fully equilibrate. The relatively short 
chain lengths considered in this study, and the relative 
low densities at which most of the phases occur, result 
in a system that does not have a pathological dynami- 
cal behavior. Therefore, the very slow pressure anneal- 
ing described above, together with the monitoring of all 
the observables considered in this study, including the 
average cluster size and its distribution where sufficient 
to establish equilibrium. Moreover, we checked that we 
could reproduce all the phases starting from completely 
different initial configurations. 

In our study we considered tethers with a minimum 
of A'^ = 5 and a maximum of A^ = 300 monomers, and 
colloids of diameter ranging from 02 ~ 2(Ti to 02 = IScti. 

All of our simulations are carried out using a total of 
512 hybrid colloids at room temperature, and the longest 
CTjS^mulations took about six months of computer time on 
an Intel Xeon X5355 2.66GHz processor. In all simula- 
tions we set the time step to At = O.OIStq. For each 
pressure annealing step, AP, simulations were run for a 
minimum of 10^ (for the small chains) to a maximum of 
10* timesteps (for the long chains.) 

Every observable reported in this paper is expressed in 
dimcnsionless units. 



III. RESULTS 

Apart from the harmonic potential, which serves a 
purely structural purpose by enforcing connectivity be- 
tween the different components of our hybrid colloid, 
there are no attractive interactions in our system. As 
a consequence, the free energy is dominated, at the low 
concentrations considered in our study, by the configu- 
rational entropy of its components. Although the con- 
figurational entropy is, strictly speaking, associated with 
hard potentials, we have chosen a large value for ks to en- 
sure that bonds are very close to their equilibrium length, 
and have also run a few simulations with a stricter ex- 
cluded volume constraint by setting ei,j = SOO/cbP. We 
find no discernible difference between the two cases under 
several thermodynamic conditions. 

What follows are the phase diagrams for ideal and self- 
avoiding tethers as a function of the volume fraction of 
colloids, (j> = Tra2Nc/ {6V), and the effective polymer- 
colloid size ratio, which we define as a = 2Rg/a2- 
Nc — 512 is the number of colloids, V the volume of 
the simulation box, and Rg is the radius of gyration of 
a polymer tether, which scales as Rg ~ (A^)^''^ for ideal 
polymers and as Rg ^ A^^/s j^j. self-avoiding ones. The 
calculated bulk prefactors for our models are 0.57 and 
0.60 for ideal and self-avoiding chains respectively. Ex- 
perimentally, one can easily control the behavior of the 
polymer by altering the properties of the solvent. For 
instance, the polymer will behave ideally at the solvent 
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9 point, and as a SAW at larger temperatures. 



Ideal chains 

Figure [2l obtained using several combinations of col- 
loidal radii and chain lengths, shows the different phases 
arising from the organization of the particles in the sys- 
tem as a function of volume fraction for different val- 
ues of a, and presents several interesting features. The 
lines are a guide to the eye, and are found by identifying 
the threshold densities and size ratios above which phase 
change occurred. 

For a sufficiently small, a < 1, the presence of the 
tethers does not alter the ability of the colloids to crys- 
tallize into a macroscopic FCC crystal once the system 
is compressed above a threshold volume fraction. This is 
exactly how tether-free colloids crystallize under analo- 
gous conditions, and is achieved in our system by chain 
localization into either the interstitial space between the 
colloids (for very small a) or into crystal vacancies as 
depicted in Fig. This is only possible as long as the 
chains are short enough to fit within a vacancy without 
exerting a significant amount of pressure arising from 
chain confinement. The formation of crystal vacancies 
is the first hint of colloidal/polymer segregation. This 
phase is preceded by a fluid phase of small micelles at 
a lower volume fraction (Fig. [5] region D). These de- 
form and freeze as the system pressure is increased into 
structurally FCC-compatible cages: the vacancies in the 
colloidal crystal lattice. Each vacancy is typically filled 
by the polymer chains of all colloids surrounding it, and 
their locations present no obvious translational order. In 
fact, we find a non-negligible number of vacancy pairs 
distributed across the colloidal crystal. 

Interestingly, for 1 < a < 1.3, the colloidal crystal 
phase ceases to form, and is replaced by a disordered 
micellar phase (see Fig. OB)- This is clearly due to the 
increased free energy cost associated with chain con- 
finement into a vacancy which grows quadratically with 
Q;.[26t Ai^ n{2Rg/ a2)^ = no? ^ where n is the number 
of chains in the same vacancy. To mitigate this effect, 
the typical cage sizes become larger and the geometrical 
rearrangement into an FCC-cell becomes expensive. The 
presence of these unstructured micelles at large volume 
fraction frustrates and disrupts the formation of a high 
density colloidal ordered phase. 

Above q; ~ 1.3, the system assembles into low-density 
micellar crystals (see Fig.[2]3, with the colloidal particles 
freely diffusing at their surfaces. This phase is analo- 
gous to that observed using a coarse-grained, soft-sphere 
model for the polymers [22|. 

We argue that the dominant contribution to the sys- 
tem pressure in the micellar regime comes from the free 
energy penalty associated with chain confinement within 
each micelle. The free energy cost per micelle associated 
with it can be readily evaluated by simple scaling argu- 



ments [1^, 

A/„„(i)\ (3) 

where is the radius of the micelle, from which we 
estimate that the pressure of the system should scale with 
micellar radius as 

P (X iVn.n^ , (4) 

where is the number of micelles forming the crystal. 

Figure [3] shows how all data collected for different com- 
binations of colloidal radii and polymer lengths in the mi- 
cellar crystal phase can indeed be collapsed into the same 
master curve when properly normalized. A power law fit 
to the data, i.e. the pressure P imposed in our simu- 
lations and the corresponding measured average micellar 
radius, yields a pressure dependance on P oc i?m*^^\ 
which is consistent with the eq. |3]for large values of i?,^. 
Clearly, our theory breaks down at very large densities, 
i.e. small micellar radii, where long tethers begin to 
radiate out of the micellar cores. This happens when 
the main mechanism of micellar shrinkage involves exclu- 
sively colloidal expulsion from the micellar surface, caus- 
ing significant thickening and layering of colloids in the 
inter-micellar regions. 

Unfortunately extending the validity of our scaling be- 
havior to large size ratios would require performing simu- 
lations with very long polymers, and this becomes quickly 
untreatable. Nevertheless, the fact that all the data gath- 
ered in our simulations for different size ratios collapse 
onto the same master curve is strongly suggestive that 
the dominant contribution to the system pressure comes 
indeed from the free energy penalty associated with chain 
confinement. 



Self-avoiding chains 

Figure H] shows the phase behavior as a function of 
particle volume fraction for different values of particle- 
to-polymer size ratio when self-avoiding chains arc con- 
nected to the colloids, and presents a quite different land- 
scape. We still find that for sufficiently small a, col- 
loids crystallize into an FCC crystal by fitting the chains 
in the colloidal interstitial spaces (Fig. UK). However, 
chains never mix to form vacancies, and as the length 
of the polymer increases, the colloidal crystal becomes 
frustrated and eventually ceases to form. Unlike the case 
of ideal polymers, we see no evidence of a micellar phase. 
We believe this is due to the large entropic barrier associ- 
ated with overlapping multiple confined chains. This can 
be estimated by computing the confinement free energy of 
a polymer of length equal to the sum of all chains in the 
cavity, which would grow as AF cx [R''^ j R^f I ^^'^ , 
where R^^ is the radius of gyration of a chain of length 
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FIG. 2: (Color online) Phase diagram of colloids with ideal tethers as a function of the polymer-colloid size ratio a and colloid 
volume fraction Snapshots of the phases in region (A), (B), and (C), depicting the colloidal crystal, the disordered micellar, 
and the micellar crystal phase, respectively, are also shown. For the sake of clarity, in snapshot (B) and (C), the colloidal 
particles are depicted using a light, low-density pixel representation, while the dark regions show where the polymer chains are 
located. 



nN and v ~ 3/5. Clearly, the free energy dependence 
on both the number of chains n and size of the cage i?„j, 
AF oc n^/^(R^/ Rm)'^-'^\ is much stronger than what ob- 
tained for ideal chains [231. ^ result, as soon as chains 
become confined, any significant amount of polymer over- 
lap is highly unfavorable. 

As the polymer size increases, for 0.5 < a < 1.75, the 
dense phase presents no colloidal order. Chains do not 
mix with each other and occupy the interstitial spaces 
in between colloids. The overall shape of the chains is 
elongated, as this geometry is entropically more favorable 
than a spherical one [28H30i |.. 

As soon as a becomes larger than 1.75, the micellar 
phase found for ideal chains is replaced by a disordered 
bicontinuous phase (Fig. HP), which allows for a more 



effective lateral packing of the chains. This phase is pre- 
ceded by the formation of small colloidal clusters driven 
together by a combination of depletion interactions and 
chain-chain repulsions (Fig. |4j3). The colloid-rich region 
presents, in both cases, a significant degree of crystalline 
order. 

The cluster phase is stable within a relatively narrow 
range of volume fractions, and is promptly transformed 
into the bicontinuous phase as soon as is sufficiently 
large for the clusters to merge. Fig. [5] shows how the 
size of the largest colloidal cluster in the system, nor- 
malized by the total number of colloids, grows with the 
system volume fraction. It is worth mentioning that col- 
loidal clusters can grow quite thick, and this can only be 
attained at the expense of the entropy of the polymers 
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FIG. 3: (Color online) Data collapse of the rescaled pressure 
P/{NmnRg) as a function of micellar size Rm for different 
values of a. 



connected to the particles at the core of each cluster, as 
they need to be partially unwrapped. We believe that 
this free energy cost may actually limit the overall thick- 
ness of the clusters and incentivize linear, rather than 
isotropic, cluster growth. 

The overall phenomenology in this region can again 
be understood in terms of chain confinement. It is well 
known [3l| that the free energy cost to completely over- 
lap two unconfined chains is about 2k^T , independent of 
the polymer length. As a result, at low volume fractions, 
there isn't a significant driving force for self-organization. 
However, as (j) increases and the chain sizes become 
smaller than i?g, IS.F acquires, as discussed above, a non- 
trivial dependence of the number of chains, n, sharing 
the same volume. This leads to chain reorganization and 
subsequent colloidal clustering. These clusters present no 
translational order or size monodispersity, and are stabi- 
lized by their mutual effective repulsions, which extend to 
a surface-to-surface range that is typically smaller than 
i?g. As soon as (j) is sufficiently large, clusters merge to 
further minimize chain-chain interactions and the bicon- 
tinuous phase discussed above is formed. For even larger 
volume fractions we observe significant ordering of the 
overall structure of the bicontinuous phase; however, the 
system sizes considered in this study are too small to 
make any conclusive claim in this regard. 

The colloidal-cluster phase can be interpreted as a dis- 
ordered inverted micellar phase. We cannot a priori ex- 
clude the existence of an inverted micellar crystal phase 
for even larger polymer lengths than the ones considered 
in this study, but such an analysis is out of the reach of 
our computational resources. 



of a colloidal particle. We find a variety of self-assembled 
structures as a function of polymer-colloid size ratio and 
volume fraction. The structures are driven by compress- 
ing the disordered low-density states and can be under- 
stood in terms of the entropy of both tethers and col- 
loids. Wc have identified chain confinement as the key 
parameter to sort out the physical mechanisms driving 
self-assembly in this system. 

It would be interesting to test whether, for self- 
avoiding polymers, an ordered bicontinuous phase and 
a crystal phase of inverted micelles can indeed be ob- 
tained, and to study how the phase behavior presented 
in this manuscript changes as a function of the number 
of grafted polymers. 

We wish to stress that both disordered and ordered 
micellar phases were observed by Capone and collabo- 
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22| while studying a system of diblock copoly- 
mers modeled as an ideal and self-avoiding polymer with 
a density-dependent effective soft-sphere potential. This 
seems to suggest that (a) the nature of the micellar phase 
for ideal tethers is not too sensitive to the details of the 
interaction, and (b) for a sufficiently large, ideal chains 
are indeed well-characterized by an additive effective pair 
potential. The problem becomes more complicated when 
dealing with self-avoiding polymers. When multiple poly- 
mers are confined within the same region their interaction 
energy does not scale linearly with the number of chains, 
but as n^/^, and up to v!^ for even larger densities [27| . 
This is clearly not pairwisc additive. Some preliminary 
results obtained using an effective soft-spherical poten- 
tial to describe the polymer (to be published elsewhere) 
indicate quite different phase behavior, including several 
ordered phases which are not found in our simulations 
with explicit polymers. This seems to suggest that a 
more sophisticated coarse-graining of self-avoiding poly- 
mers is required to obtain the correct phenomenological 
behavior of this system. 
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IV. CONCLUSIONS AND DISCUSSION 



We report the phase behavior of a system of hybrid col- 
loids formed by grafting a single polymer on the surface 
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FIG. 4: (Color online) Phase behavior for the case of self-avoiding chains. In the graph in the top-left corner of the figure, the 
vertical axis indicates the polymer-colloid size ratio q and the horizontal axis is the colloidal volume fraction (j). Snapshots of the 
phases in region (A), (B) and (C), depicting the colloidal crystal, the colloidal cluster, and the bicontinuous phase respectively 
are also shown. For the sake of clarity, in snapshots (B) and (C), the polymers are depicted using a light, low density pixel 
representation. 
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